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ABSTRACT 

The mobility of atoms, molecules and radicals in icy grain mantles regulate ice restructuring, des¬ 
orption, and chemistry in astrophysical environments. Interstellar ices are dominated by H 2 0, and 
diffusion on external and internal (pore) surfaces of H 2 O-rich ices is therefore a key process to con¬ 
strain. This study aims to quantify the diffusion kinetics and barrier of the abundant ice constituent 
CO into H 2 0 dominated ices at low temperatures (15-23 K), by measuring the mixing rate of ini¬ 
tially layered H 2 0(:C0 2 )/C0 ices. The mixed fraction of CO as a function of time is determined by 
monitoring the shape of the infrared CO stretching band. Mixing is observed at all investigated tem¬ 
peratures on minute time scales, and can be ascribed to CO diffusion in H 2 0 ice pores. The diffusion 
coefficient and final mixed fraction depend on ice temperature, porosity, thickness and composition. 

The experiments are analyzed by applying Fick’s diffusion equation under the assumption that mixing 
is due to CO diffusion into an immobile H 2 0 ice. The extracted energy barrier for CO diffusion into 
amorphous H 2 0 ice is ^160 K. This is effectively a surface diffusion barrier. The derived barrier is 
low compared to current surface diffusion barriers in use in astrochemical models. Its adoption may 
significantly change the expected timescales for different ice processes in interstellar environments. 

Keywords: astrochemistry; molecular processes; methods: laboratory; ISM: molecules 


1. INTRODUCTION 

In the cold (T<20 K) regions of the interstellar 
medium, the surfaces of dust grains are coated by ice 
due to a combination of freeze-out of gas-phase molecules 
and an active grain surface chemistry. Based on obser¬ 
vations of ice absorption bands in protostellar and cloud 
lines-of-sight, the main ice constituent is H 2 0 followed 
by CO and C0 2 (see lGibb et al.ll2Q04l: lOberg et al.ll 201 lL 
for reviews). In most lines of sight there is spectroscopic 
evidence for separate H 2 Q- and CO-rich ice ph ases, with 
C0 2 mixed into both (jPontoppidan et aLll2008l ). The dif¬ 
fusion efficiency of molecules between these two phases is 
poorly constrained, limiting our understanding for how 
the ice morphology evolves when the grains are heated 
during star formation. For example, the extent of dif- 
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fusion of volatile species into H 2 0 ice will regulate the 
importance of entrapment of volatiles in H 2 0 ice. 

Diffusion in ices is also important for the chemi¬ 
cal evolution during star formation. Ices are major 
reservoirs of volatiles and ice chemistry is the pro¬ 
posed pathway to both simple volatiles such as H 2 0 
and CH 3 OH, and to the complex orga nic m olecules ob¬ 
served toward some protostars (e.g. Tielens & Hagenl 
119821 : IGarrod et al.ll 20 Q 8 l : iHerbst fc van Dishoeckl 120091 ). 
The efficiencies of these pathways crucially depend on 
the mobility of th e reactants on ice surfaces and in the 
bulk of the ices (IGarrod fe Paulvl 120111 : iGarrodl 120131 : 
IVasvunin fe Herbstl 120131 ). 

Up until recently there were few laboratory con¬ 
straints on the diffu sion of molecules on top of or inside 
of in t erstellar ices ([L ivingston et al. 2002 ; Smith et all 
120101 : iMispelaer et al.1 120131 : iKarssemeiier et al. 2 014 ). 
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Most astrochemical models therefore parameterize sur¬ 
face and bulk diffusion barriers as fractions (30-80%) 
of the b etter understood mole c ular desorption barri¬ 
ers (e.g., iTielens fc Hagenl 119821: iGarrod fe Paulvl [20111 : 
iChang fc Her bst 2012]). Surface diffusion is generally 
modeled as a hopping process between different poten¬ 
tial mini ma with a specific barrier or ensemble of barri¬ 
ers (e.g. JChang et all200 5: lCuppen et al.lll2013l ). Within 
the ice matrix, pore walls provide internal surfaces and 
diffusion in these pores can be treated similarly to ex¬ 
ternal surface diffusion. ‘Proper’ bulk diffusion has been 
imagined as either a swapping p rocess I 


_ ^ __ _ Derg et al. _ 

iFuchs et "all 120091: lGarrodll2013[) . or a m ovement into in¬ 
terstitial spaces between molecules (e.g.. lLamberts et~ahl 
I2013L12014 iChang fc Her bstlfe 0141 ). Both processes have 
been modeled to have high barriers compared to surface 
diffusion. 

CO diffusion in H2O ice has been the focus of sev¬ 
eral studies, both because of its interstellar relevance 
and because its utility as a model system when eval¬ 
uating different experimental approaches. lOberg et al.l 
(200 9) constrained the segregation rate of thin (10s of 
monolayers) pre-mixed hGOiCO ices at 23-27 K, un¬ 
der ultra-high vacuum conditions, and obtained a bar¬ 
rier for CO diffusion out of H2O-CO ice mixt ures of 
300±100 K (26±9 meV). iMispelaer et al.l (20131 ) found 
a diffusion barrier of 120T170 K (10T15 meV) for CO 
diffusion through a thick (100s of monolayers) H2O ice 
film a nd into the gas-phase. In a similar set of experi¬ 
ments. iKarssemeiier et al.l (20141 ) measured the diffusion 
of CO out of COiIUO ice mixture, through a thick amor¬ 
phous ice layer, and into the gas-phase at temperatures 
of 32-50 K and found a diffusion barrier of 300T170 K 
(26=bl5 meV). The experimentally determined barrier of 
CO diffusion out of H 2 0-rich ices of 120-300 K (10- 
26 meV) is somewhat lower than the >400 K (34 meV) 
CO diffusion barrier in current use in astrochemical mod¬ 
els (Garrodll2013l ). 

The experimental results are in some tension with re¬ 
cent s imulations of CO d i ffusio n on amorphous H2O, 
where IKarssemeiier et ahl (2 014 ) found that there are 
two populations of CO binding sites, one described as 
a strong-binding nanopore site and the other as a weaker 
surface site, with diffusion barriers of 80 meV and 30 
meV, respectively. Experiments and simulations may be 
reconcilable if diffusion kinetics in H2O ice mainly de¬ 
pends on the barrier height of the weakest bound CO. To 
isolate the low-barrier diffusion from diffusion between 
more strongly bound sites, we study the diffusion kinet¬ 
ics of CO in H2O ice at very low temperatures, where 
CO in strongly bound sites will be completely immobile. 
The study of low-temperature diffusion is facilitated by 
our experimental strategy focusing on CO diffusion into 
H2O ice rather than out of the ice. The latter must 
be performed at temperatures above the CO desorption 
temperature of ^30 K. 

In summary we aim to quantify the diffusion rate of 
CO by measuring the mixing rates of initially layered 
H2O/CO ices. §2 presents the new experimental set-up 
designed to explore the physics and chemistry of thin 
ices, and the experimental procedures and spectral anal¬ 
yses specific for this study. The modeling strategies are 
explained and motivated in §3. §4 presents the measured 


mixing kinetics of initially layered H2O/CO ices and their 
dependencies on ice temperature, thickness, morphology 
and composition. The experimental results are modeled 
using a realization of Fick’s diffusion equation to extract 
the CO diffusion barrier. The results are discussed in 
light of previous experimental and theoretical work, and 
astrochemical modeling in §5 and summarized in §6. 

2. EXPERIMENTAL DETAILS 
2.1. Experimental set-up 



Figure 1 . Three-dimensional view of the chamber, with the main 
experimental instrumentation marked. 


The experiments were all carried out in a new labo¬ 
ratory set-up (Fig. [I]), designed for interstellar ice ana¬ 
log experiments. It consists of a 13” (330.2 mm) spher¬ 
ical stainless steel ultrahigh vacuum (UHV) chamber 
(custom-made, Pfeiffer Vacuum), evacuated by a Pfeif¬ 
fer Turbo HiPace 400 pump backed by a DUO 10M ro¬ 
tary vane pump to a base pressure of ^ 10“ 10 mbar at 
room temperature. The ices are grown on a 2 mm thick 
IR transparent Csl substrate with a 19 mm clear view, 
mounted on an optical ring sample holder at the center of 
the chamber. The sample holder is connected to the cold 
tip of a closed cycle He cryostat (Model CS204B, Ad¬ 
vanced Research Systems, Inc.) capable of cooling the 
Csl substrate down to 11 K. The cryostat is mounted on 
the top port of the chamber via a differentially pumped 
UHV rotary seal (Thermionics RNN-400) that allows 360 
degree rotation of the Csl substrate inside the cham¬ 
ber without breaking the vacuum. The Csl substrate 
is mounted onto the nickel plated OHFC copper sample 
holder using silver gaskets for good thermal contact. A 50 
ohm thermofoil heater is installed on the cryocooler tip so 
that the temperature of the substrate can be varied be¬ 
tween 12-350 K. The substrate temperature is regulated 
by a temperature controller (LakeShore Model 335) using 
two calibrated silicon diode sensors (accuracy of 0.1 K), 
one connected directly to the sample holder and the other 
near the heater element. 

Ice composition, thickness and morphology are mon¬ 
itored through infrared absorption transmission spec¬ 
troscopy at 4000-400 cm -1 using a Bruker Vertex 70v 
spectrometer with a liquid nitrogen cooled MCT detec¬ 
tor. The IR beam from the spectrometer enters and exits 
the UHV chamber through KBr ports located on either 
side of the chamber. The IR beam is focused onto the 
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Csl substrate by a custom designed set of transfer optics 
(flat and off-axis paraboloidal mirrors) located inside the 
interferometer chamber, and refocused onto the MCT de¬ 
tector using a similar optics assembly located inside the 
detector chamber. The spectrometer, the external inter¬ 
ferometer chamber and the MCT detector chamber are 
evacuated to 2 mbar to avoid atmospheric interference 
with the ice spectra. The spectra reported in this paper 
have a resolution of 1 cm -1 and have been background 
subtracted. A Pfeiffer quadrupole mass spectrometer 
(QMG 220M1, mass range 1-100 amu, and resolution of 
0.5 amu) is positioned 40 mm off the Csl substrate. It is 
used to continuously monitor the gas composition during 
ice growth and warm-up of the deposited ices. 

Ices are grown in situ by exposing the cold Csl sub¬ 
strate to a constant flow of gas from an independently 
pumped gas-line (base pressure lower than 10 -5 mbar) 
using a gas doser consisting of an xyz stage (enabling 
positioning 1-90 mm from the substrate window), a 
precision leak valve (MDC vacuum), a deposition tube 
(4.8 mm diameter, 14 inches (355.6 mm) long). The de¬ 
position tube is connected to a gas-mixing line, which 
is differentially pumped down to 7 x 10 -4 mbar, moni¬ 
tored with baratron and pirani gauges. The line contains 
various ports for glass bulb vessels or stainless steel can¬ 
isters, which can be interchanged as needed, and several 
on/off valves, one of which allows the gases to flow to the 
precision leak valve on the gas doser. 

2.2. Experimental procedure 

The experiments were carried out with CO (>99%, 
GT & S, Inc.), C0 2 (>99%, Isotec), deionized water 
and deuterated water (99.9%, Cambridge Isotope Lab¬ 
oratory, Inc.). The water and deuterated water samples 
were purified by three freeze-thaw cycles under vacuum. 
Pure gas samples were prepared by filling one of the bulbs 
on the gas-mixing line to a pressure of 10 mbar. 10 mbar 
gas mixtures were made by adding species sequentially 
to one of the bulbs at the correct proportions, estimated 
using the pressure gauges and and the known volumes 
of the different bulbs and tubes within the gas mixing 
line. This procedure results in maximum ~10% devi¬ 
ations from target mixtures, based on QMS measure¬ 
ments. 

Prior to each experiment the gas flow was set to 
5 x 10 -7 mbar at room temperature. The Csl window 
was then cooled down to 12 K, followed by the intro¬ 
duction of a sequence gases or gas mixtures to build up 
a layered ice structure. The deposition tube was kept 
close to the surface to minimize porosity in the case of 
deposition along the surface normal (most experiments) 
or a highly controlled porosity structure for the handful 
of ices deposited at an angle. 

Figure [2] depicts a general summary of the various ex¬ 
perimental sample preparations. To obtain a layered ice 
structure, CO was deposited first with an ice thickness 
of 5-60 monolayers (ML). The ice growth was controlled 
using a known gas flow rate (measured at room temper¬ 
ature) and the deposition time. The achieved CO ice 
thickness was then measured by acquiring an infrared 
spectra normal to the sample of the deposited ice us¬ 
ing the relationship between column density, integrated 
optical depth and band strength: 


where Ni is the column density (molec/cm 2 ), and 
f Ti{iy)du is the integrated optical depth of the IR band 
area, and A{ is the band strength of the species i. The 
CO ice thickness in ML is calculated using the extracted 
CO column density and the known CO ice density. The 
CO ice layer is then covered by a layer of H 2 0 ice (or a 
H 2 0 ice mixture), whose thickness is regulated and mea¬ 
sured using the same procedure. Table Q] presents the 
band strengths, band integration range and ice densities 
of the species used in the different experiments. 


Table 1 

Ice infrared and density data. 


Species 

IR band 

Int. range 
cm -1 

Ai 

cm molec -1 

P 

g cm -3 

CO 

C-O str. 

2120-2170 

l.lxl0 -17a 

0.81 b 

h 2 o 

O-H str. 

3000-3600 

2.0x 10 _16a 

0.94 c 

H 2 0:C0 2 

O-H str. 

3000-3600 

1.6xl0 _16d 

0.94 e 

co 2 

C-O str. 

2310-2370 

7.6xl0 _17a 

1.3 e 

d 2 o 

O-D str. 

2225-2700 

1.4xl0 16f 

1.04 c 

^Gerakines et al. (1995|), fc LoeffIer et al. (20051. 

iJenniskens & Blake (1994), d Oberg et al. (2007|), 

lEscribano et al. (2013), IVenvaminov & Prendergast (1997) 


Following the initial ice characterization the sample 
was quickly heated (5 K/min) to the targeted experi¬ 
ment temperature and was then maintained at this tem¬ 
perature (15-23 K) for 1.5-4.5 hours while monitoring 
changes in the ice morphology using the infrared spec¬ 
tral features of CO and H 2 0. The first IR spectrum was 
collected immediately upon reaching the target temper¬ 
ature (t = 0) and subsequent IR spectra were collected 
every two minutes with 64 scans per spectrum. Dur¬ 
ing the entire experiment the gas-phase composition was 
also monitored using the QMS, checking especially for 
ice desorption - no CO desorption was observed during 
any of the experiments. CO desorption was also checked 
by monitoring the total CO IR band area during each 
experiment. We find that the total band area is con¬ 
served and thus that CO desorption is negligible also at 
the highest experiment temperature of 23 K. 

2.3. IR spectral band analysis 

Pure CO ice and CO mixed with H 2 Q have distinct 
spectral bands (e.g.. iSandford fe Allamandolall19881 Fig. 
[3^, b). Pure CO exhibits a Lorentzian IR band profile 
centered at 2139 cm -1 ([Bouwman et al.ll2007lh while the 
spectral feature of CO in a CO:H 2 0 ice mixture can be 
fit by two overlapping Gaussians, centered at 2137 cm -1 
and 2152 cm -1 . The 2152 cm -1 band has been iden¬ 
tified with CO molecules interacting with dangling OH 
sites within the water ice, and is therefore often referred 
to as the ‘polar’ band. The origin of the 2137 cm -1 
band is less well understood. It has been attributed 
to CO interactions with water molecules whose hydro¬ 
gen atoms are bound within the ice matrix, and also 
to CO interacti ons with other CO molecules in H 2 Q 
ice environment ([Schmitt et al.l 119891: iManca et al.ll2Q0ll 
lAl-Halabi et al.ll2004D . 
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Figure 2. An outline of the different ice configurations used on the presented experiments. The top arrows indicate the deposition angle; 
d indicates the interface of CO and water, and h indicates the vacuum/ice film interface. In most experiments a thin layer of CO ice was 
covered by H 2 O ice (a). In a sub-set of experiments the CO layer was instead covered by a H20:C02 ice mixture (b). In a third set of 
experiments the CO was covered by a layer of H 2 O followed by a third layer of D 2 O. All these ices were deposited along the surface normal 
to minimize porosity. A small set of ices were instead deposited at an angle to increase the ice porosity (d). In each case we then monitored 
the mixing of CO into the ice over layers. 



2170 2160 2150 2140 2130 2170 2160 2150 2140 2130 2120 

Wavenumbers / cm 1 


Figure 3. Infrared spectra of pure CO ice together with a fitted Lorentzian profile (a), a H^OiCO 5:1 ice mixture (b) together with the 
spectral Gaussian fits, the changing spectral features as a H 2 O/CO layered ice mixes (c), and the fit of pure and mixed CO ice components 
to the end result of (c) using the Lorentzian and Gaussian functions from (a) and (b). 


Figure [3b shows the evolving spectra of an initially 
layered H 2 0/C0 ice kept at 17 K for 90 min. The ini¬ 
tial spectra can be fit as a superposition of a pure CO 
ice spectra and a CO ice mixture spectra, with the pure 
CO ice dominating. With time the spectra approach the 
profiles of a completely mixed ice. The final spectrum 
(Fig. [3ji) is well fitted by two Gaussians, and a small 
Lorentzian. This 3-component fit was automated using 
the IDL MPFIT function and applied to all experiments 
and times. Because of the frequency overlap between 
the pure CO band and the 2137 cm -1 mixed ice band, 
the area of the 2137 cm -1 band is less well constrained 
compared to the 2152 cm -1 band. Since the 2152 cm -1 
band correlates with mixed CO in water ices, we devel¬ 
oped an analysis method using this band without having 
to isolate the pure band from the 2137 cm -1 mixed ice 
band. To validate this approach we explored the relative 
band intensities of the 2137 cm -1 and 2152 cm -1 bands 


in five separate C0:H 2 0 ice mixtures using our exper¬ 
imental set-up. The mixtures spanned mixing ratios of 
1:2 to 1:10 and also included one experiment with C0 2 
mixed in. At 12 K all CO mixture spectra present a 1- 
to-3 intensity ratio between the 2152 cm -1 feature and 
the total CO band (Fig. 0]). The 2152 cm -1 band alone 
thus provides a good measure of the total amount of CO 
mixed with H 2 0. We quantify the amount of mixed CO 
at each time step in the diffusion experiments with, 


A" mix 


^2152 X 3 

^total 


( 2 ) 


where 7V m i x is the mixed fraction of CO, A 2 i 5 2 is the 
integrated area of the 2152 cm -1 feature, and A tot ai is 
the total integrated area of the CO stretch band. 


3. DIFFUSION MODELING 
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Figure 4. The integrated intensity ratio of the total CO band 
area over the 2152 cm -1 (polar) spectral feature in different ice 
mixtures. 


d h 



Figure 5. Theoretical progression of CO-H 2 O through CO diffu¬ 
sion when starting with layered ice films of CO and H 2 O. In this 
example h— 10 nm, d— 1 nm, and D = 10 -16 cm 2 s _1 . 


We use Fick’s second law of diffusion to model the 
observed C0:H 2 0 ice mixing as a function of time. In a 
one dimensional system, Fick’s second law describes the 
concentration of a diffusing species, c(z,£), as a function 
of time, £, and position, z: 


dc(z, t ) 
dt 


D(T) 


d 2 c(z , t ) 
dz 2 


(3) 


This law should apply to the C0:H 2 0 mixing kinet¬ 
ics if the mixing is dominated by the diffusion of the 
more volatile CO into the H 2 0 matrix described by the 
temperature dependent diffusion constant D(T). The 
equation wa s previously applied to C O diffusion ex¬ 
per iments by iKarssemeijer et ahl ( 20141 ); iMispelaer et, al. i 1 
ml . where CO was modeled as diffusing through a 
H 2 0 matrix followed by desorption into the vacuum. 

In the context of Fick’s diffusion equation the initially 
layered H 2 0/C0 ice system is described such that the 
substrate is at z = 0 , z = d marks the interface between 
CO and the H 2 0 layer, and the surface of the H 2 0 layer 
is at z = h. In this study, the boundary conditions for 
solving Fick’s equation is set by the lack of CO desorp¬ 
tion, which entails that the total amount of CO in the 
system is constant and therefore, the flux of CO at the 
substrate (z = 0) and vacuum (z = h ) interfaces is 0, i.e., 

dC Q Z z ^ = 0 at z = 0 and z = h. In addition the CO con¬ 
centration was initially assumed to be Co between z = 0 
and z = d, and zero elsewhere. These initial and bound¬ 
ary conditions give the following solution to Equation [3l 



The generic solution as a function of z and t is shown 
in Fig. [5l As function of time, CO spreads out into the 
H 2 0 matrix, approaching the final mixed state, i.e., at 
t = 0 , all CO is located from 0 to d, and with time more 
and more CO diffuses across the interface d, approaching 
complete mixing at infinite times. 

The mixed fraction of CO, A mix , is represented as a 
function of time by integrating Eq. 0]over the H 2 0 layer 
and dividing by the total amount of CO: 


1 f h 

Wnix(£) = / c (z, t)dz 

dco Jd 


h — d 
h 


2 h . 0 (nnd\ 

n= 1 x 7 


x exp 


h 2 



( 5 ) 


Based on Fig. [3t, the initial stage of the mixing exper¬ 
iments is not well-described by a complete lack of mix¬ 
ing, however. To allow for mixing on deposition of the 
H 2 0 ice layer, and mixing during the fast warm-up from 
the deposition temperature to the mixing temperature at 
t < 0 a time offset to was added to Eq. O We also added 
a nuisance parameter No to account for experimental un¬ 
certainties in the measured ice thicknesses, yielding 


Nmi x(t) 


N 0 


h — d 
h 


E 


2 Noh . 2 
■ sm 


n 2 7r 2 d 



x exp 


h 2 


-D(t + to) 


When using this equation to fit experiments, h and d 
are taken from spectroscopic measurements, and D, No, 
and to are modeled as free parameters. Fig. [ 6 ] shows the 
application of this solution to layered H 2 0/CO experi¬ 
ments of different thicknesses (Exps. 2 and 13 in Table 
El, demonstrating how the H 2 0/CO mixing is modeled 
as a function of time including a measured amount of ice 
mixing on deposition. The model solution can be fit to 
H 2 0/CO ices with any thickness and mixing ratio, but 
when the H 2 0 ice is thin compared to the CO:H 2 0 in¬ 
terface, there is no well-defined H 2 0 matrix for the CO 
to diffuse into, and the derived diffusion coefficients are 
difficult to interpret. That is, the extracted diffusion co¬ 
efficients will no longer be applicable for diffusion in pure 
H 2 0 ice. 


4. RESULTS 

4.1. H 2 O/CO ice mixing dependencies 

Table [2] lists the experimental details for all initially 
layered H 2 0/CO ice mixing experiments. The exper¬ 
iments cover a range of ice mixing temperature ( 12 - 
23 K), CO and H 2 0 ice thicknesses (5-59 ML and 13- 
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Table 2 

H 2 O/CO ice mixing experiments ordered based on the main parameter varied within each subset of experiments, together with initial and 
final fit values and the diffusion rate from Fick’s diffusion law analysis in §3.3. Reported uncertainties do not include systematic ones (see 

text). 


Exp. 

H 2 0/CO 

(ML) 

C^dep 

n 

Tmix 

(K) X 

JV mix (0) b 

(%) 

^mix(^final) 

(%) 

D 

(cm 2 s -1 ) 

Ice temperature 







1 

34/10 

0 

12 

26 

26 

3.8x 10 -18 

± 5.8 x 10 -18 

2 

31/8 

0 

15 

49 

64 

5.2x 10 -17 

± 1.4 x 10 -18 

3 

35/9 

0 

16 

51 

65 

3.0X10 -16 

± 1.1 x 10 -17 

4 

28/8 

0 

17 

44 

67 

3.1 x 10 -16 

±4.6 x 10 -18 

5 

31/9 

0 

20 

40 

64 

2.1 x 10 -15 

± 6.5 x 10- 17 

6 

29/10 

0 

23 

37 

61 

2.1 x 10 -15 

± 5.1 x 10 -17 

Ice thickness 







7 

13/8 

0 

15 

46 

55 

2.1 x 10 -16 

± 1.3 x 10 -17 

8 

59/10 

0 

15 

38 

63 

1.4 x 10 -16 

± 5.4 x 10 -19 

9 

14/9 

0 

20 

53 

56 

9.2 x 10 -15 

± 1.3 x 10 -15 

10 

59/9 

0 

20 

44 

88 

2.0 x 10 -15 

± 3.4 x 10- 16 

11 

37/5 

0 

17 

68 

88 

2.6x 10 -16 

±4.7 x 10- 17 

12 

40/14 

0 

17 

11 

42 

6.1 x 10 -16 

± 1.3 x 10 -18 

d 2 o/h 2 O/CO 







13 

89/81/14 

0 

15 

16 

56 

1.3x 10 -16 

± 6.1 x 10 -18 

14 

117/106/30 

0 

15 

8 

43 

1.8x 10 -16 

± 2.9 x 10 -18 

15 

107/86/59 

0 

15 

5 

23 

7.7x 10 -16 

± 1.6 x 10 -17 

Porosity a 







16 

30/9 

30 

15 

49 

70 

8.5x 10 -16 

± 8.9 x 10- 17 

17 

14/9 

45 

15 

49 

66 

6.4x 10 -16 

± 1.4 x 10 -16 

18 

19/8 

30 

20 

45 

68 

3.3x 10 -15 

± 3.9 x 10 -16 

19 

26/8 

45 

20 

56 

81 

4.0x 10 -15 

± 9.6 x 10 -18 

Reproducibility 







20 

34/10 

0 

17 

26 

56 

2.2xl0 -16 

±4.7 x 10 -18 

21 

39/10 

0 

17 

30 

61 

2.9x 10 -16 

±4.6 x 10 -18 

22 

35/10 

0 

17 

32 

55 

2.6x 10 -16 

± 2.1 x 10 -18 

H 2 0:C0 2 /C0 







23 d 

46/10 

0 

15 

28 

48 

2.0xl0 -16 

± 5.7 x 10 -18 

24 d 

34/9 

0 

17 

29 

57 

1.5x 10 -16 

± 2.1 x 10 -18 

25 d 

35/10 

0 

20 

24 

68 

6.0x 10 -16 

± 5.2 x 10 -18 


a Porosity of the water layer was varied by changing the angle of incidence, «dep? during deposition of the the H 2 O. 
b Presented in both the percentage of CO that had mixed and the corresponding monolayers that had mixed. The uncertainty for each of these 

values was <10%. 

c Mixing rate, fcdiff? was determined from D divided by the square of the water (top) layer thickness. 
d The top layer was a mixture of 4:1 H 2 0:C02- 



Height z (nm) 

Figure 6. The modeled concentration of CO with respect to ice 
height (z), at various time stamps, for thin and thick ice experi¬ 
ments. 

223 ML, respectively), H 2 O/CO abundance ratios (1.5- 
14), H 2 O ice porosities (using different deposition an¬ 
gles), and compositions of the top H 2 O-rich ice layer. 
The fiducial experiment with a H 2 O/CO ratio of 4/1 
and a mixing temperature of IT K was repeated four 
times during the experimental series to test the experi¬ 
mental reproducibility. 

Figure 0 shows the time progression of CO mixing with 


the H 2 O ice for a representative subset of the H 2 O/CO 
experiments listed in Table [2j Above 12 K all layered 
ice experiments begin to mix within the first few min¬ 
utes, and the final mixed fraction is typically reached 
within tens of minutes. In Fig. [7^, the mixing rate, 
most clearly seen from the initial slope of the curves, in¬ 
creases with temperature, while the final mixed fraction 
does not. In ice experiments with different ice thickness 
(Fig. [7)3— d), the initial mixed CO fraction depends, as 
expected, on the thickness of the CO under layer. The 
final mixed CO fraction increases with the thickness of 
the H 2 O layer when the CO underlayer thickness is held 
constant, and decreases with the CO underlayer thick¬ 
ness when the thickness of the H 2 O ice layer is constant, 
i.e. it seems to increase with increasing H 2 O/CO ratio. 
At 15 K, the mixing rate is higher for ices deposited at 
an angle compared to ices deposited along the surface 
normal (black vs. blue curves in Fig. [Tt), while at 20 K, 
the segregation is too fast to visually infer any difference 
in mixing rate between ices deposite d at different a ngles . 
Ba sed on pr evio us experiments (e.g. JKimmel et ahlhOOli : 
iRaut et aP 120071 ). an increasing deposition angle results 
in a more porous H 2 O ice when the ice is deposited at 
low temperatures. The experiments thus show that the 
CO mixing rate increases with ice porosity. 

To quantify these dependencies we fit all experiments 
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Figure 7. The fraction of mixed CO ice vs. time in initially layered CO/H 2 O ices. Panel a) shows the increasing mixing rates with 
increasing mixing temperature, Exps. 1, 2, 4 & 5. Panel b) and c) shows the mixing dependencies on H 2 O ice thickness covering 10 ML 
of CO at 15 (Exps. 7, 2 &; 8) and 20 K (Exps. 9, 5 & 10), respectively. Panel d) shows the mixing dependency on CO thickness at 17 K 
when covered by 35-40 ML of H 2 O ice, Exps. 11, 4 & 12. Panel e) and f) show CO mixing with water layers deposited at different angles 
away from the surface normal; Exps. 2, 16 &; 17 at 15 K, and Exps. 5, 18 &; 19 at 20 K. 


using our solution to Fick’s diffusion equation (Eq. [6j). 
The resulting fits are shown in Fig. [8j and the extracted 
initial and final ice mixing data, and the diffusion coef¬ 
ficients are reported in Table [2] The fits to the data are 
generally excellent for the experiments, resulting in small 
fit uncertainties. As outlined below the experimental+fit 
uncertainties and the uncertainties in experiments char¬ 
acterized by fast mixing are considerably larger. 

There is always a significant amount of mixing at the 
beginning of each experiment. Between 5 and 68% of 
the CO starts out mixed with the H 2 0 (dependent on 
the CO ice thickness), corresponding to a pre-mixed in¬ 
terface 2-4 ML thick. The mixed CO fraction at the end 
of each experiment, A%i x (£fi na i), varies between 26 and 
88%, and increases with increasing H 2 0/CO ratio as ex¬ 
pected from theory: the expected final mixing ratio at 
t = oc is No(h — d)/h. This is best exemplified when 
comparing experiment 9 and 10, and 11 and 12, and 14 
and 15. 

The extracted diffusion coefficients depend strongly on 
temperature between 12 and 20 K; it increases by three 
orders of magnitude over this temperature range for the 
H 2 0/CO experiments with a thickness of ~32/9 (above 
23 K, the diffusion is too fast to effectively measure with 
our time resolution). For the same ice thickness the 
fft+experimental uncertainty (lcr) is estimated to 20% 
based on four repeated experiments at 17 K. 

The measured diffusion coefficients are expected to be 
independent off ice thickness as long as diffusion into pure 
H 2 0 ice dominates and the mixing time scales are long 
compared to the time resolution. The former condition 
depends on a combination of ice thickness and H 2 0/CO 
ratio. The mixing time scale decreases with increasing 
temperature and with decreasing ice thickness. The thin 
ice experiments at high temperatures thus result in the 
highest model uncertainties. In experiments with short 
mixing time scales, the extracted diffusion coefficients are 
lower limits (especially experiments 6-7, 9-10, 18-19). 


Empirically, measured diffusion coefficients at 15, 17 and 
20 K do not depend on ice thickness (within 50%) as long 
as the H 2 0/C0>3 and the total ice thickness is >36 ML; 
cf exps. 2, 8, 13 and 14, and exps. 4 and 11. The 
thickest ices in these comparisons contains a third D 2 0 
layer on top of the H 2 0 layer and it is assumed that this 
isotopic substitution does not affect the CO diffusion into 
the ice. As inferred visually, the diffusion coefficients do 
depend on ice porosity. At 15 K, the diffusion coefficient 
for the ice deposited at an angle (exp. 16) is an order of 
magnitude higher compared to a the diffusion coefficient 
derived for a comparable ice deposited along the surface 
normal (exps. 2 and 8). 

In addition to the pure ice layer experiments we car¬ 
ried out a few experiments, inspired by interstellar ice 
compositions, with C0 2 mixed into the H 2 0 layer. In 
these experiments, ^40 ML of a 1:4 C0 2 :H 2 0 mixture 
was layered on top of ~10 ML of CO, and the tempera¬ 
ture was varied from 15 to 20 K. While the initial mixed 
fraction was constant for all three experiments, both the 
diffusion coefficient and the final mixed fraction increased 
with temperature, with the coefficient tripling between 
15 to 20 K and the final mixed fraction increasing from 
0.48 to 0.68. (Fig. [9]). The increase in the diffusion coef¬ 
ficient with temperature in these experiments is smaller 
compared to the H 2 0/C0 mixing experiments without 
C0 2 as a layer component. 

4.2. CO diffusion kinetics and barriers 

We use the diffusion coefficients derived in §4.1 to con¬ 
strain the diffusion barrier assuming an Arrhenius type 
process. For a single barrier process D(T) should follow 
the Arrhenius equation, 

D = T x e~ Editt/T , (7) 

where T is a pre-exponential factor, the diffusion 
barrier in K, and T the ice temperature. If D is known for 
multiple ice mixing experiments performed at different 
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Figure 8. Solutions of the Fickian model (solid blue lines) for all experiments (red circles). 



Figure 9. The mixed CO fraction versus time in the CC^IEO ice 
mixture experiments (Exps. 23-25). A 17 K H 2 O/CO experiment 
(Exp. 4) is shown for comparison. 

temperatures, E^iff is readily extracted by fitting a linear 
function to ln(Z}) versus 1/T. The pre-exponential fac¬ 
tor can also be extracted from these fits, but its physical 
meaning is difficult to interpret without additional con¬ 
straints on the mixing process. Figure [T0h shows that the 
diffusion coefficients from Exps. 1-6, 20-22 (where all 
parameters except for temperature were kept constant) 
are well-fit by the Arrhenius equation. The resulting ice 
mixing barrier is 158T12 K. The reported uncertainties 
in D in Table [2] are generally very small because sys¬ 
tematic uncertainty is not incorporated. Based on fits 
to repeated experiments the experimental uncertainty is 
closer to ~20%, and when using the derived diffusion 


coefficients to calculate barriers, a 20% uncertainty was 
therefore added to each data point. 

FigurefTOb.c shows the fit of the Arrhenius equation to 
other subsets of experiments where temperature is the 
only variable; i.e., ‘thin’ (H 2 O/CO 14/8 ML) and ‘thick’ 
(H 2 O/CO 59/10 ML) experiments, more porous ices de¬ 
fined by the deposition angle, and the experiments with 
H20:C02 ice mixture layered on top. Because of the 
small number of data points in each case the fit param¬ 
eters are highly uncertain, but the CO 2 mixture exper¬ 
iments and possibly the high porosity experiments do 
have significantly lower mixing barriers compared to the 
other experiments (Table [3j). The pre-exponential factors 
are also reported. 


Table 3 

Calculated energy barriers and pre exponential factors assuming 
an Arrhenius-type mixing process. 


Experiment 

r 

(cm 2 sec -1 ) 

Ediff 

(K) 

Ediff 

(meV) 

Mixing temperature 

3.1 Xl0 -12 

158±12 

14 

Thin water layer 

7.9 xlO -10 

227±26 

20 

Thick water layer 

5.8 XlO -12 

159±25 

14 

Porosity, 30° 

1.9 XlO -13 

81±27 

7 

Porosity, 45° 

9.9 xlO -13 

110±28 

9 

CO 2 mixture 

1.7 xlO -14 

72±18 

10 
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Figure 10. Arrhenius plots incorporating the diffusion coefficients 
from Exps. 1-6 and 20-22 are plotted in panel a), the diffusion 
coefficients from the water thickness (Exps. 7 & 9, 8 & 10) and 
porosity (Exps. 16 & 18, 17 & 19) experiments plots in b) and c), 
and the H20:C02 ice mixtures (Exps. 23-25) in panel d). 


4.3. The CO diffusion mechanism: constraints from 
H 2 O ice spectroscopy 

The H 2 0 ice (and its isotopologu e) infrare d band pro¬ 
file changes in the pre sence of CO (jRowland et al.l 1991; 
iBonwman et al.l I2007D , which can provide independent 
constraints on the ice mixing rate and mechanism. These 
changes are most pronounced on the blue wing of the 
H 2 0 stretching band, where OH-bonds that are not in¬ 
teracting with the hydrogen bonding network in the bulk 
of the water ice absorb (Fig. fTTh). In pure H 2 0 ice, these 
dangling OH bonds trace pore surface area, with a double 
peak IR profile at 3719 and 3697 cm -1 . H 2 0 molecules 
that interact with CO present a similar feature slightly 
shifted to the red (jBouwman et al.ll2007[ ). at 3639 cm -1 . 
There are corresponding dangling OD bands in pure and 
mixed D 2 0 ices, with the double pure features at 2750 
and 2731 cm -1 , and the mixed feature at 2692 cm -1 . 

Figure fTTh and b shows that during the diffusion ex¬ 
periments there is a clear loss of the H 2 0 3719 and 
3697 cm -1 features and a simultaneous growth of the 
3639 cm -1 band. Figure ITTl l shows the corresponding 
time series, which displays a similar profile to the mixing 
time series based on CO spectroscopy, but with a lower 
signal-to-noise ratio. This increasing exchange of H 2 0 
pore spectral features with C0-H 2 0 spectral features in¬ 
dicates that CO diffuses into the H 2 0 ice through pores. 

To constrain the diffusion mechanism further we run 
several experiments with a D 2 0 layer on top of the H 2 0 
layer and simultaneously monitored the changes in the 
H 2 0 and D 2 0 dangling OH and OD bands. Figure ITTll 
and e shows that when the bottom CO layer is thin 
(14 ML), no CO diffusion into the top D 2 0 layer is ob¬ 
served. As the CO ice layer thickness is increased, CO 
diffusion into the D 2 0 layer is observed, but always at 
a lower level compare to the intermediary H 2 0 layer. 
The initial CO mixing rate also appears lower for D 2 0 
compared to H 2 0 as would be expected since CO has 
to travel through H 2 0 before entering into the D 2 0 ice. 
The low SNR in these experiments as well as an un¬ 
known amount of D/H exchange reactions between D 2 0 
and H 2 0 prohibits an independent measure of the diffu¬ 
sion rate based on these experiments. Figure fl2l shows, 
however, that the data qualitatively compares well with 
the Fickian model solutions for these experiments, based 
on the fit to the CO spectral time progression modeling 
n §4.2.2. CO diffusion into a H 2 0 matrix as described 
by our solution to Fick’s diffusion law thus seems to be 
a good model of the H 2 0/C0 ice mixing process. 

5. DISCUSSION 

The presented experiments aimed at characterizing CO 
(surface) diffusion in H 2 O-rich ices by measuring the mix¬ 
ing of initially layered ices and extracting the diffusion 
coefficient from the data. This approach should work as 
long as the observed ice mixing is dominated by CO dif¬ 
fusion into the H 2 0 ice matrix. The general agreement 
between data and models based on Fick’s diffusion law 
supports that CO diffusion indeed drives the mixing pro¬ 
cess. Without additional information the CO diffusion 
could be due either bulk diffusion though e.g., molecular 
swapping (jOberg et al.ll2009D. or hopping on the s urfaces 
of nano- and micropores (jKarssemeiler et ahll2014l ). Both 
could a priori result in the observed diffusion profiles. It 
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Figure 11. Spectra of pure H 2 O and a H 2 0:C0 mixture displaying the differences in the dangling OH/OD spectral regions (panel a). 
Panel b and c show the time evolution of the dangling OH/OD difference spectra in an ice mixing experiment, and panel d and e show the 
integrated intensity increase with time in the three D 2 O/H 2 O/CO experiments (Exps. 13-15). 

is furthermore not necessary that the same diffusion pro¬ 
cess dominates at all times, since mixing of CO (or other 
molecules) into the ice could change both the hopping 
and swapping barriers. 

There are several experimental and theoretical consid¬ 
erations that suggest that diffusion of CO in our ice sys¬ 
tems is dominated by pore diffusion, however. First, the 
derived diffusion barrier is very low (<170 K). Second, 
the diffusion rate increases by an order of magnitude at 
15 K when the ice porosity is increased. This may be 
partially due to a reduction in the diffusion barrier (see 
below), but a more straightforward explanation is that 
most of the diffusion rate increase is due to an increase in 
pore area, and that the diffusion process in these exper¬ 
iments proceeds similarly to the low-porosity ice experi¬ 
ments. Finally and perhaps most importantly, the H 2 O 
ice spectroscopy in these experiments reveals that as the 
CO diffuses into the H 2 O ice, the dangling OH bonds is 
replaced by spectroscopic features assigned to dangling 
OH-CO interactions, which is expected if the pore walls 
become increasingly covered by CO. 

Simulations provide additional clues on the CO-H 2 O 
in teractions i n these pores. In particular, simulations by 
iKarssemeiier et ahl ( 2014 ) show that water ice surfaces 
contain multiple sites that interact strongly or weakly 
with CO. These interactions should manifest themselves 
in different CO spectral features in CO:H 2 0 ice mix¬ 
tures, but exact identifications are still lacking. At our 
low experimental temperatures, CO that becomes bound 
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Figure 12. The Fickian model of CO diffusion into H 2 O and D 2 O 
when the D 2 O is layered on top. The model reproduces both the 
delay into the D 2 O layer compared to the H 2 O layer, and how the 
final portioning of CO between the H 2 O and D 2 O layers approaches 
unity as the CO under layer thickness is increased. 
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in a strongly interacting site is expected to be become 
bound and the diffusion kinetics should be regulated 
by the number density and the diffusion barrier of the 
weakly interacting sites. As the strongly bound sites 
fill up, the CO can diffuse further into the ice before 
it gets trapped, consistent with the diffusion behavior in 
the thick D 2 O/H 2 O experiments. Spectroscopically, we 
measure mixing by the growth of CO interacting with 
dangling OH, which probably correspond to CO that is 
trapped in a strongly bound site, or nanopore. In this 
scenario the derived diffusion coefficient thus describes 
the rate at which a molecule diffuses through the ice 
along weakly bound sites, while the remaining parame¬ 
ters in the diffusion equation characterizes the availabil¬ 
ity of empty nanopores/strongly bound sites to diffuse 
into. In this scenario, the diffusion length scale of a CO 
molecule depend on the number of weakly bound sites be¬ 
fore a strongly bound site is encountered. The traveled 
distance should thus increase as more sites become occu¬ 
pied by CO molecules. The maximum diffusion length is 
limited by the amount of available CO. This is confirmed 
by the isotopically layered thick ice experiments where it 
is clear that is is shown that if the CO under layer is 
thin, no CO will make it all the way through the H 2 0 
ice into the D 2 0 ice. Figure [13] shows a schematic of 
the proposed mixing mechanism of the CO and H 2 0 ices 
through CO diffusion into the H 2 0 ice via micropores 
and nanopores. 

The experiments further provide evidence for that the 
already low CO diffusion barrier of ^160 K is further re¬ 
duced if the ice porosity is increased and if C0 2 is mixed 
into the H 2 0 ice. The simplest explanation for the low 
diffusion barrier in the high-porosity experiments is that 
CO has access to a higher number of low-barrier diffu¬ 
sion pathways. It may also be due to a hi g her density 
of ‘la rge’ micropores (jKimmel et al.l 120011: iRaut et ah! 
Hqq 3), i-e., larger than a few molecules across, which 
could enable CO molecules to “piggy-back” on top of 
other CO molecules coating the micropore walls. In other 
words, the ^100 K barriers derived in the high-porosity 
ice experiments may reflect the CO-CO diffusion barrier. 
More experiments are required to distinguish between 
these scenarios, and indeed to ensure that the extracted 
small barrier is real. The C02:H 2 0 ice mixture exper¬ 
iments both result in a lower barrier and in an overall 
diffusion rate. This suggests that the addition C0 2 both 
changes the overall binding environment, adding lower 
binding sites compared to pure H 2 0 ice, and reduces the 
number of pores that CO can diffuse into. 

5.1. Astrophysical implications 

Based on the presented experiments and analysis, the 
barrier for CO to diffuse into H 2 0 ice along pore walls 
is <170 K, with a best fit value of ^158 K for a pure, 
low-porosity ice. This compares well, within the uncer¬ 
tainties, with results from previous studies on CO dif¬ 
fusion in amorphous H 2 0 ice using three different ap¬ 
proaches: segregation of CO from CO:H20 ice mixtures 
(jOberg et al.1 l2009|h CO des orption from CO:H 2 0 ice 
mixtures (|Karssemeijer et al.ll2014l). and CO des orption 
from layered H 2 O/CO ices (jMispelaer et al.ir2013[ ). Sum¬ 
marized in Tabled] the derived CO diffusion barriers in 
these studies all range between 120 and 300 K (10-30 


meV). Considering the good agreement between these 
very different experiments (including large differences in 
vacuum and deposition conditions, ice thicknesses and ice 
morphology), it appears robust that there are CO diffu¬ 
sion barriers that are significantly lower than currently 
assumed in astrochemical models for either surface or 
bulk diffusion. 

Table 4 

Summary of our and previous CO diffusion results. 


Diffusion Barrier 
K meV 


Experiments 


Our experiments 

158±12 

Mil 

Segregation 

300±100 

26±9 

Out of thick H 2 O layer b 

116=bl74 

10±15 

Out of CO:H20 mixture 0 

302=L174 

26=bl5 

Models 



Gas-Grain d 

400 

34 

CO free on H 2 O surface 0 

348 

30 

CO trapped in H 2 O surface pore° 

929 

80 


IQberg et al.l (120091) . ^Mispelaer et al.l (120131) . 
3Kar ssemeiier et al. (120141) . ^Garrod ( 2 0131) 


Our measured CO diffusion barrier describes CO diffu¬ 
sion in H 2 0 ice with pores, i.e. CO diffusion on internal 
H 2 0 ice surfaces. In non-porous ice, the derived barrier 
cannot be used to describe ice bulk diffusion, only dif¬ 
fusion on the ice surface. Furthermore, unless the H 2 0 
is impure, resulting in some of the strongly bound sites 
being occupied by other volatiles, the diffusion length 
scale with this low barrier may not be very long. Most 
astrophysical ices may fall into this category based on 
the form ation process o f the H 2 0 ice in the interstellar 
medium (Ga rrodl 120131) and its observed intimat e mix¬ 
ing with other volatiles iPontoppidan et al.l (12008 1. The 
derived barrier thus provides an important constraint on 
CO surface diffusion; the H 2 0-CO diffusion barrier of 
^160 K is a factor of 2-3 lower than what is currently 
assumed in astrochemical models for surface diffusion. 
This discrepancy implies that diffusion of CO, and per¬ 
haps many other species, is possible at much lower tem¬ 
peratures than currently assumed. 

6. CONCLUSIONS 

The kinetics of CO diffusion into water ice at low tem¬ 
peratures (below the CO desorption temperature) were 
examined using initially layered ices. Based on the diffu¬ 
sion dependencies on temperature, ice thickness, poros¬ 
ity, and composition and the subsequent analysis we con¬ 
clude that: 

1. CO is mobile in low-porosity amorphous H 2 0 ice 
at 15-23 K, with a diffusion length scale that de¬ 
pends on the number of available strongly bound 
sites where CO can become trapped. 

2. The measured CO diffusion into the H 2 0 ice ma¬ 
trix is a pore-mediated process and characterizes 
CO diffusion on H 2 0 surfaces, rather than bulk 
diffusion. 

3. The observed dependence of the CO diffusion on 
ice temperature implies that the lowest CO diffu- 
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c) 


Figure 13. Schematic of the mixing process. The layered system at t = 0 is shown in a). Represented in b), mixing occurs, with the 
inset showing the CO molecules diffusing along the micropore surfaces into the strong-binding nanopore sites. After some period of time, 
the layer becomes fully mixed, shown in c). 


sion barrier into the pores of amorphous H 2 O ice 
is 160 K . 

4. This barrier is lower in more porous ices and when 
adding CO 2 is added to the H 2 O matrix. 

5. The derived barrier is low compared to existing val¬ 
ues in astrochemical networks, indicative of that 
surface processes are more efficient than currently 
assumed. 


We gratefully acknowledge productive discussion with 
E. Herbst, D. Pesce, S. Liss, and P. King. D. Graninger 
assisted with several of the experiments. This work was 
supported by OT2 funding from Herschel Space Observa¬ 
tory (PI I. Cleeves), European Research Council (ERC- 
2010-StG, Grant Agreement no. 259510-KISMOL), 
and the VIDI research program 700.10.427, financed by 
The Netherlands Organization for Scientific Research 
(NWO). 

REFERENCES 

Al-Halabi, A. and Fraser, H. J. and Kroes, G. J. and van 
Dishoeck, E. F. 2004, A&A, 422, 777 
Bouwman, J., Ludwig, W., Awad, Z., et al. 2007, A&A, 476, 995 
Chang, Q., Cuppen, H. M., & Herbst, E. 2005, A&A, 434, 599 
Chang, Q. & Herbst, E. 2012, ApJ, 759, 147 
Chang, Q. & Herbst, E. 2014, ApJ, 787, 135 
Cuppen, H. M., Karssemeijer, L. J., Sz Lamberts, T. 2013, 
Chemical Reviews, 113, 8840 
Escribano, R. M., Munoz Caro, G. M., Cruz-Diaz, G. A., 
Rodriguez-Lazcano, Y., & Mate, B. 2013, Proceedings of the 
National Academy of Science, 110, 12899 
Fuchs, G. W., Cuppen, H. M., Ioppolo, S., et al. 2009, A&A, 505, 
629 

Garrod, R. T. 2013, ApJ, 765, 60 
Garrod, R. T. & Pauly, T. 2011, ApJ, 735, 15 
Garrod, R. T., Weaver, S. L. W., & Herbst, E. 2008, ApJ, 682, 
283 


Gerakines, P. A., Schutte, W. A., Greenberg, J. M., & van 
Dishoeck, E. F. 1995, A&A, 296, 810 
Gibb, E. L., Whittet, D. C. B., Boogert, A. C. A., Sz Tielens, 

A. G. G. M. 2004, ApJS, 151, 35 
Herbst, E. Sz van Dishoeck, E. F. 2009, ARA&A, 47, 427 
Jenniskens, P. & Blake, D. F. 1994, Science, 265, 753 
Karssemeijer, L. J., Ioppolo, S., van Hemert, M. C., et al. 2014, 
ApJ, 781, 16 

Kimmel, G. A., Stevenson, K. P., Dohnalek, Z., Smith, R. S., Sz 
Kay, B. D. 2001, J. Chem. Phys., 114, 5284 
Lamberts, T., Cuppen, H. M., Ioppolo, S., Sz Linnartz, H. 2013, 
Phys. Chem. Chem. Phys., 15, 8287 
Lamberts, T., de Vries, X., Sz Cuppen, H. M. 2014, Faraday 
Discussions, 168, 327 

Livingston, F., Smith, J., Sz George, S. 2002, J. Phys. Chem. A, 
106, 6309 

Loeffler, M. J., Baratta, G. A., Palumbo, M. E., Strazzulla, G., Sz 
Baragiola, R. A. 2005, A&A, 435, 587 
Manca, C., Martin, C., Allouche, A., Sz Roubin, P. 2001, J. Phys. 
Chem. B, 105, 12861 

Mispelaer, F., Theule, P., Aouididi, H., et al. 2013, A&A, 555, 
AT3 

Oberg, K. I., Boogert, A. C. A., Pontoppidan, K. M., et al. 2011, 
ApJ, 740, 109 

Oberg, K. I., Fayolle, E. C., Cuppen, H. M., van Dishoeck, E. F., 
Sz Linnartz, H. 2009, A&A, 505, 183 
Oberg, K. I., Fraser, H. J., Boogert, A. C. A., et al. 2007, A&A, 
462, 1187 

Pontoppidan, K. M., Boogert, A. C. A., Fraser, H. J., et al. 2008, 
ApJ, 678, 1005 

Raut, U., Fama, M., Teolis, B. D., Sz Baragiola, R. A. 2007, 

J. Chem. Phys., 127, 4713 

Rowland, B., Fisher, M., & Devlin, J. P. 1991, J. Phys. Chem., 

95, 1378 

Sandford, S. A. Sz Allamandola, L. J. 1988, Icarus, 76, 201 
Schmitt, B., Greenberg, J. M., Sz Grim, R. J. A. 1989, ApJ, 340, 
L33 

Smith, R. S., Matthiesen, J., Sz Kay, B. D. 2010, J. Chem. Phys., 
133, 174504 

Tielens, A. G. G. M. & Hagen, W. 1982, A&A, 114, 245 

Vasyunin, A. I. Sz Herbst, E. 2013, ApJ, 762, 86 

Venyaminov, S. Sz Prendergast, F. 1997, Ana. Biochem., 248, 234 





























